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A SEMIGRAPHICAL METHOD FOR A DEORBIT TRAJECTORY DESIGN STUDY 


By Richard N. Green, William F. Hampshire H, 
and Sue W. Souders 
Langley Research Center 

SUMMARY 

A semigraphical method to aid in the selection of a deorbit trajectory which is com- 
patible with arbitrarily selected mission constraints is presented. The method is appli- 
cable to the case where the initial orbit about the planet and the deorbit trajectory are 
coplanar. Basically, the method consists of examining a family of candidate trajectories 
and eliminating, by a graphical cross-plotting technique, those which violate the mission 
constraints. In addition, the trade-offs which exist between the various demands on the 1 
trajectory become graphically obvious. The method is flexible in that the user is free to 
dictate the mission mode and the mission constraints. An example illustrating the use of 
the method for a Mars mission is presented. 

INTRODUCTION 

Man has long sought to understand and explore the solar system around him. In 
recent years, considerable knowledge of the neighboring planets has been obtained from 
interplanetary probes. These missions have consisted of instrumented space vehicles 
which passed close to the planets while proceeding in orbits about the Sun. Therefore, 
experiments to obtain planetary information were short in duration. A likely candidate 
for future exploration involves placing satellites in orbits about the planets whereupon 
experiments can be performed over long periods of time. Surface probes can be ejected 
from such a satellite orbit to gather more comprehensive planetary data. This procedure 
requires a velocity maneuver to transfer the probe to a descent trajectory. Therefore 
the combined use of the orbiting satellite and the ejected surface probe seems a likely 
candidate for use in interplanetary exploration in the near future. The descending probe 
may be required to satisfy a number of constraints imposed by the mission goals, the 
subsystem requirements, or by the laws of physics. In the final result, these considera- 
tions must be integrated into a workable and mutually compatible spacecraft system. 
Eventually, the spacecraft orbit and the probe trajectory must be designed to satisfy the 
mission objectives within the constraints imposed by the entire set of subsystems. There 
is a requirement, therefore, to determine a method of examining the many facets of such 
a complex mission. 



The intent of this report is to present a semigraphical method to aid in the selection 
of descent trajectories which are compatible with the mission constraints. To reduce the 
complexity of the general three-dimensional problem, the motion of the probe has been 
assumed to remain in the plane of the satellite orbit. This assumption would seem to be 
reasonable since considerably more fuel would have to be expended to establish an out -of - 
plane trajectory. Basically, the method consists of examining a family of reasonable tra- 
jectories and eliminating, by a graphical cross -plotting technique, those which violate the 
mission constraints. In addition, the trade-offs which exist between the various demands 
on the trajectory become graphically obvious. A discussion of the method and its use for 
various mission modes is given. Finally, the method is applied to a hypothetical Mars 
mission and results in the identification of a family of descent trajectories which satisfy 
all the imposed constraints over a range bounded by two Martian atmospheric models. 

SYMBOLS 

A cross-sectional area, meters2 

Cj} drag coefficient 

f true anomaly of deorbit, degrees {see fig. 1(a)) 

h altitude of probe, meters {see fig. 1(c)) 

h a apoapsis altitude of satellite orbit, kilometers (see fig. 1(a)) 

altitude of probe at parachute deployment, meters 
hp periapsis altitude of satellite orbit, kilometers (see fig, 1(a)) 

l television footprint length, meters (see fig. 1(c)) 

m mass, kilograms 

m/Cp)A ballistic coefficient of probe before parachute deployment, kilograms/meter^ 

(m/CpjA)^ ballistic coefficient of probe -parachute system after parachute deployment, 
kilogram s/meter2 

V velocity of probe, meter s/second (see fig. 1(b)) 
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AV magnitude of deorbit velocity maneuver, meter s/second (see fig. 1(a)) 

|3^ angle at probe measured from local vertical to line of sight to satellite, 

degrees (see fig. 1(b)) 

y flight-path angle of probe, degrees (see fig. 1(b)) 

>' y of probe at limiting altitude of atmosphere, degrees 

Ay, maximum change in y, during parachute descent phase, degrees 

6 application angle of AV measured from satellite negative velocity vector 

at deorbit, degrees (see fig. 1(a)) 

angle at probe measured from probe's negative velocity vector to line of 
sight to satellite, degrees (see fig. 1(b)) 

y^ ^ y^ at parachute deployment, degrees 

v l min minimum y^ during parachute descent phase, degrees 

$ position angle of probe at impact referenced to periapsis of satellite orbit 

and considering no atmosphere during descent, degrees 

position angle of probe at impact referenced to periapsis of satellite orbit 
and considering an atmosphere during descent, degrees (see fig. 1(a)) 

ip any relative position angle between probe and satellite 

\p^ lower constraint on \p 

1^2 upper constraint on ip 

SEMIGRAPHICAL METHOD 

For the purpose of illustrating the method, consider the following mission mode and 
and mission constraints. The mission mode consists of a velocity maneuver out of orbit, 
aerodynamic deceleration, and parachute braking near the terminal phase. The parachute 
is to be deployed at a specific Mach number. To insure a survivable impact on the 
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surface of the planet, the altitude at which deployment occurs is restricted to be above 
some minimum altitude h^. Assume that the communication system requires that a rel- 
ative position angle ^ between the probe and the satellite be bounded at impact for the 
purposes of transmitting surface data to the satellite. Therefore, the position angle at 
impact is constrained between and These constraints could cover a number of 

different situations and are defined here in general terms for the purpose of illustration. 
The size and shape of the orbit have been defined along with the position in orbit at which 
the deorbit maneuver is to be performed. In addition, the aerodynamics of the probe and 
the parachute system have been defined. The problem reduces to bounding the magnitude 
and direction of the velocity maneuver to produce a family of impacting trajectories 
which will permit the successful accomplishment of the mission and will simultaneously 
satisfy, or not violate, any of the imposed mission constraints. 

Since the probe remains in the plane of the satellite orbit, the descent trajectory 
can be uniquely defined by two independent parameters once the initial orbit and the deor- 
bit point have been fixed. These two independent parameters are of the utmost impor- 
tance in the graphical representation of the deorbit geometry. The easiest ones to 
visualize are AV, the magnitude of the deorbit velocity maneuver, and 8 , the angle of 
application (fig. 1(a)), since these are the control variables during the actual mission. 
Nevertheless, several different sets of parameters were examined, and by trial and error 
another set was found to be more useful for design purposes, namely, 6 and the vacuum 
landing point on the surface. Consider sketch (a) where $ is the vacuum landing point on 
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f, h a , and h p fixed 

AV 


>e 


the surface measured from some reference line, taken here as the line of apsides of the 
initial orbit, and y e is the flight-path angle at atmospheric entry. Lines of constant AV 
and y e have been plotted as a function of $ and 8 . Notice that f, the true anomaly at 
deorbit; h a , the apoapsis altitude; and h p , the periapsis altitude have been fixed. Since 
the trajectory can be defined by <t> and 8 , the sketch is a graphical representation of 
all possible descent trajectories. It remains to eliminate those areas which violate the 
mission constraints. These areas were eliminated by selecting a number of sample 
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trajectories that encompassed ail possibilities. These sample trajectories corresponded 
to the intersections of the lines of constant AV and y e . For each sample trajectory 
the true landing point was computed by numerically integrating the equations of 
motion through the atmosphere. The resulting landing point is displayed in sketch (b) 

f, h a , and hp fixed 
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Sketch (b) 

where the symbols represent the data points. Several of the sample trajectories have the 
same flight-path angle at atmospheric entry whereas other groupings have common values 
of AV at deorbit. This result is due to choosing the samples at the intersections of the 
constant y e and AV lines in sketch (a). By connecting the appropriate data points, 
lines of constant AV and y e were formed on the new coordinate plane. Thus , the 
AV,y e grid was plotted as a function of 8 and an integrated parameter The 

approach of plotting the AV,y e grid against different parameters is the basis of the 
method presented here. In this manner the true landing point of all possible trajectories 
was found by examining only a few well chosen samples. 

Consider the mission constraints to eliminate the undesirable trajectories. For 
example, the parachute system is restricted to deployment above h<j to insure a surviv- 
able impact. The computed altitude of parachute deployment for each sample trajectory 
was then plotted as a function of 8 and the AV,y e grid was constructed (sketch (c)). 
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f, h a , and h D fixed 



Sketch (c) 
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The trajectories within the shaded area violate the parachute deployment constraint and 
were eliminated from consideration. In a similar manner the position angle at impact 
was plotted as a function of (sketch (d)). Only those trajectories which have a if; 
at impact between pi and ^2 are acceptable from a communication point of view. 



Sketch (d) 

All constraint lines were then transferred to sketch (e) by means of the AV,y e 
grid which was common to all plots. In this way a family of descent trajectories was 


f, h a , and hp fixed 
* = *1 
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Sketch (e) 

found which satisfied the mission objectives. Within this area, additional optimization can 
be made; that is, the most desirable descent trajectory can be selected. If the mission 
is over constrained and all trajectories are eliminated, one must reexamine the various 
constraints to find a solution. This examination is immediately possible with the sepa- 
rate parametric plots (sketches (c) and (d)) without generating additional trajectory data. 




To reiterate, a family of desirable descent trajectories was found for a given mis- 
sion mode. The vacuum landing point and the deorb it velocity application angle were 
taken as the independent parameters to define the trajectory, A group of sample trajec- 
tories were chosen to represent various combinations of AV and y e and were then 
numerically integrated through the atmosphere. The undesirable trajectories were elimi- 
nated by examining the separate parametric plots associated with each constraint. The 
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constraint lines were then transferred to a common plot by means of the AV,y e grid. 

By this semigraphical method a family of descent trajectories which satisfied the mission 
constraints was found. 


DISCUSSION OF METHOD 

Normally, the mission mode and constraints are more complex than those in the 
illustrative example. Although only two constraints were considered, any number of con- 
straints could have been imposed. Provided a parameter is single-valued along each tra- 
jectory, one can plot the parameter as a function of some appropriate secondary parameter 
and the AV,y e grid constructed. The secondary parameter should be chosen on the basis 
of the shape of the resulting AV,y e grid. It has been found that this choice has a great 
effect on the graphical representation of the data. The mission mode could have been 
extended to include a more sophisticated deployment criteria for the parachute. The use 
of thrusting rockets to slow the probe could have been included in the study. In any event, 
the procedure should prove to be applicable for more complex missions. 

Many different parameters can be investigated by repetition of the method. The true 
anomaly which was held constant in the previous example is a case in point. The influence 
of the deorbit position on the selection of a descent trajectory could be investigated by 
repeating the same procedure with a different true anomaly. The family of satisfactory 
descent trajectories which result from the second iteration could be overlaid with those 
from the first. The result of this type of procedure would possibly lead to a smaller fam- 
ily of trajectories which satisfies the mission constraints over a region of deorbit true 
anomalies instead of only at one point, A similar study could be performed to investigate 
the influence of the size and shape of the satellite orbit. If the planet's atmosphere is not 
well defined, one may use the method to select a descent trajectory which insures the mis- 
sion objectives over a range of atmospheric models. 

Obviously, the method is very flexible in that the user is free to define the mission 
mode and the various constraints imposed on the descent trajectory. In addition, the 
method can be repeated to investigate the influence of various parameters. The difficulty 
in applying the method arises in the formulation of a given constraint so that it can be 
handled graphically as outlined. Many times this formulation will require some insight 
into the problem and some additional computations. More complex mission constraints 
are examined in the section devoted to the selection of a Mars descent trajectory. Because 
of the flexibility of the method, it should provide a logical, orderly procedure for a wide 
variety of studies. 
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NUMERICAL EXAMPLE 


To illustrate the method further, a study to design the descent trajectory for a hypo- 
thetical Mars-probe mission is presented. The satellite and the probe are injected into 
an elliptical orbit about Mars from which the probe will separate from the satellite and at 
an appropriate time transfer to a descent trajectory by means of a small coplanar velocity 
maneuver. The probe is slowed down by using both atmospheric braking and a parachute 
to provide for a survivable impact. While the probe is descending on the parachute, a 
television system gathers imagery data of the surface and transmits it back to the satellite. 
Constraints pertaining to both the communication and the television subsystems are con- 
sidered. In addition, the descent trajectory is chosen to insure a successful mission over 
a range of possible atmospheric models. The following numerical values are postulated 
for use throughout the numerical study, 

(a) Periapsis altitude of the orbit, hp, 1000 km 

(b) Apoapsis altitude of the orbit, h a , 20 000 km 

(c) True anomaly of deorbit, f, 240° 

(d) Voyager Mars atmospheric models, VM-8 and VM-4 

(e) Limiting altitude of the atmosphere, 243 840 m (800 000 ft) 

(f) Ballistic coefficient of the probe before parachute deployment, m/CjjA, 

39kg/m2 (0.25 slug/ft2) 

(g) Ballistic coefficient of the probe -parachute system after parachute deployment, 

(m/C D A) d , 3.8 kg/m2 (0.024 slug/ft2) 

(h) Parachute is deployed at a Mach number of 1.6 

(i) Parachute is released at an altitude of 1524 m (5000 ft) 

(j) Television lens full field of view, 32° 

(k) Television coverage of the planet begins when l = 9144 m (30 000 ft) 

(see fig, 1(c)) 

(l) Antenna on probe has a beam width of 60° 

The following mission constraints are imposed on the descent trajectory: 

(a) -20° S Y e = -15° 

(b) AV £ 400 m/s 

(c) o° ^ e ^ 180 ° 

(d) Altitude at which television coverage is initiated g 6096 m (20 000 ft) 
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(e) -30° £ i>£ £ 30° during parachute descent phase 

(f) j3j at impact £ -60° 

(g) /3^ at impact plus 5 minutes £ 60° 

The flight-path angle at atmospheric entry is bounded to avoid possible "skip- out" 
or excessive heating. The weight allocation has restricted the amount of fuel for the 
deorbit maneuver. The application angle of the velocity maneuver 9 is applied so that 
the descent trajectory lies inside of the spacecraft orbit. If the velocity increment were 
applied in an outward direction, the probe would pass outside of the satellite orbit and then 
Would cross the orbit before landing on the surface. This condition is avoided to eliminate 
possible communication problems. To avoid oblique imagery, the television coverage is 
initiated at an adequate altitude. In order to transmit these pictures to the satellite with 
a directional antenna, the descent trajectory is chosen so that the satellite will be in the 
proper position relative to the lander during the parachute descent. A measure of the 
satellite's position relative to the probe's body axis is the i> ^ angle (fig.l(b)). Since the 
probe's antenna beam width is 60°, the angle during the parachute descent is bounded 
between ±30°. Low data rate communication is insured during the entire descent phase 
and for 5 minutes after impact. To insure that the communication is not broken because 
of the horizon or possible mountain ranges, the angle (fig. 1(b)) at impact is greater 
than -60°. Similarly, the (3^ angle at 5 minutes after impact is not to exceed 60°. 

These mission objectives and mission constraints are considered in the selection of a 
descent trajectory which will provide a successful mission within the range of Martian 
atmospheres assumed to be bounded by the VM-4 and VM-8 atmospheres. 

The vacuum plot (fig. 2) was constructed to select the sample trajectories. At this 
point the constraints on y Q and 0 were considered, that is, sample trajectories were 
chosen so that y e was between -15° and -20° while 9 was bounded between 0° and 180°. 
In addition, each sample had a AV of 400 m/s or less. Therefore, a group of sample 
trajectories was chosen which encompassed all reasonable candidates. Notice that the 
vacuum plot is independent of the atmospheric model and the mission mode because of the 
geometric and conic relationships between the four parameters, 9, y , and AV. 

For this reason, the sample trajectories were applicable for the consideration of both the 
VM-8 and VM-4 atmospheres. The thinner VM-8 atmosphere was considered first. The 
sample trajectories were numerically integrated through the Martian atmospheres to 
obtain the parameters associated with each mission constraint. The true landing point 
which is dependent on both the atmospheric model and the mission constants such as the 
ballistic coefficients is presented in figure 3 (a). The altitude constraint for the initiation 
of the television system was investigated with figure 3(b). Those trajectories which vio- 
lated the two angle constraints became obvious with the construction of figures 3(c) 
and 3(d). The constraint on the relative positions of the probe and satellite for 
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communication purposes was satisfied by consideration of figure 3(e), This plot will be 
discussed subsequently. The mission constraints as listed previously were satisfied by 
transferring all the constraint lines to a common plot and by eliminating those areas of the 
trajectories which violated any of the constraints. Therefore, the shaded area in fig- 
ure 3(f) represents the family of descent trajectories that satisfied all the mission require- 
ments in a VM-8 atmosphere. 

The v i angle constraint restricted the descent trajectory to be chosen so that the 
v £ angle was maintained between ±30° during the parachute descent phase. This con- 
straint was not handled as the others were because the parameter of interest was bounded 
during a period of time and not at one specific point. Thus, a different approach was taken. 
Examination of the v ^ angle showed that the angle started decreasing immediately after 
parachute deployment and continued decreasing until a minimum was reached at which 
time it started increasing slowly until the end of the parachute phase was reached. The 
maximum value always occurred at deployment whereas the minimum value occurred at 
some other time during the parachute deceleration phase. Several different graphical 
representations were tried before figure 3(e) was considered. The abscissa is the value 
of at parachute deployment and the ordinate is the maximum variation of during 
the parachute phase. The vertical line at ^ ^ = 30° is obviously one extreme of the 
usable area of trajectories. The other can be reasoned as follows. If the v i angle at 
deployment is 20°, then the total change in vj must be less than 50° (Ay, = -50°) 
to insure that the minimum is greater than -30°. This approach establishes the 
constraint line labeled " v ^ m ^ n = -30°." These two lines bounded the area of descent tra- 
jectories which satisfied the constraint that -30° =»i = 30° during the parachute phase. 
From a communication viewpoint, the preferable condition would be equal angular excur- 
sions. This condition exists along the line labeled "v^ balanced'” ^°^ ce that this con- 
straint was examined on a plot that contained neither <5^ nor 8 as the others did. 
Nevertheless, the constraint line was transferred to a common plot by use of the AV,y e 
grid. 

The family of descent trajectories that satisfied all the mission requirements in the 
VM-4 atmosphere is represented by the shaded area in figure 4. The procedure followed 
was the same as previously outlined with the exception of the use of a different atmospheric 
model. The two areas from figures 3(f) and 4 were then transferred to the vacuum plot 
(fig. 5) since it was independent of the atmosphere. The common area, therefore, repre- 
sents the descent trajectories that satisfy the mission constraints and provide for a suc- 
cessful mission over a range of possible Martian atmospheres bounded by VM-4 and VM-8 
atmospheres. 

One of the major advantages of the semigraphical method is the presentation of the 
various trade-offs that exist. Within the common area of descent trajectories, a choice of 
a specific trajectory can be made. If one desires a steep flight -path angle, the magnitude 


10 



of the deorbit velocity is increased and the probe is forced to land farther from the peri- 
apsis of the satellite orbit. IE very shallow entry angles y e are desired the /3j con- 
straint at impact plus 5 minutes could be relaxed. This trade-off between flight-path 
angle and communication time after impact and other similar trade-offs are byproducts of 
the constraint plots. These comparisons can be seen with no additional computation or 
plotting. It is very easy to over constrain the descent trajectory; in which case, these 
trade-offs become increasingly important. 

CONCLUDING REMARKS 

A semigraphieal method to aid in the selection of a descent trajectory has been pre- 
sented. The method is flexible in that the user dictates both the mission mode and the 
mission constraints, A wide variety of studies can be performed by repeating the method 
for different mission parameters such as the dimensions of the satellite orbit, the position 
within the orbit of the deorbit maneuver, or the atmospheric model. The method expedites 
the selection of a descent trajectory since an orderly procedure has been established and 
a descriptive set of design plots have been found. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 19, 1968, 

789-30-01-01-23, 
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(a> Orbital geometry. 
Figure 1.- Mission geometry. 







(a! Variation of true landing point with deorbit A V angle* 

Figure 3,- Mars deorbit trajectory envelope for a probe descent through the VM-8 atmosphere* h p = IGOO Km; b a - 20 000 Km : f = 24$; m/CpA - 39 kg/m 2 10*25 slug/ft 2 ) 
(m/CoA) d = 3J7 kg/m 2 10*024 slug/ft 2 ); parachute deployment Mach number, 1.6; parachute release altitude, 1524 meters 15000 ft)* 






s, deg 


{b\ Variation of altitude of probe at l = 9144 meters DO ODD ft) with deorbit AV angle. 
Figure 3.- Continued. 
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(c> Variation of ^ with at impact 

Floure 3 - Continued. 
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Figure 4.- Mars deorbit trajectory envelope for a probe descent through the VfYM atmosphere. h p = 1000 km : h a = 20 000 km; f = 24CP; m/C^A - 39 kg/m? 10.25 slug/ft?) 
(m/CpA) j = 3.77 kg/m? fO,Q24 slug/ft2) ; parachute deployment Mach number, 1.6, parachute release altitude, 1524 meter's <5000 ft). 
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